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Abstract 


The main focus of this thesis is to present a numerical study of Fredholm integral equations 
of the nonlinear integro-differential type. This includes examining both the regular and 
weakly singular cases, as well as the fractional case. To obtain numerical solutions for these 
equations, we use popular projection methods like the Galerkin method and the collocation 
method, along with classical orthogonal polynomials. The primary benefit of this approach 
is that it allows us to transform the main equations for each case into a nonlinear algebraic 
system. We can then use iterative methods to solve these systems efficiently. To show the 
accuracy and effectiveness of our approach, we present several numerical examples throughout 
the thesis. These examples demonstrate how our numerical process accurately solve the given 
equations, which further confirms the effectiveness of our proposed method. 

Keywords: Orthogonal polynomials. Fredholm integro-differential equations. Galerkim 


method. Collocation technique. Nonlinear equation 


Resumé 


L’objectif principal de cette thése est de présenter une étude numérique des équations in- 
tégrales de Fredholm de type intégro-différentiel non linéaire. Cela comprend l’examen 
a la fois des cas réguliers et faiblement singuliers, ainsi que du cas fractionné. Pour 
obtenir des solutions numériques pour ces équations, nous utilisons des méthodes de 
projection populaires comme la méthode de Galerkin et la méthode de collocation, ainsi 
que des polynémes orthogonaux classiques. Le principal avantage de cette approche est 
qu’elle nous permet de transformer les équations principales de chaque cas en un sys- 
téme algébrique non linéaire. Nous pouvons ensuite utiliser des méthodes itératives 
pour résoudre efficacement ces systémes. Pour démontrer la justesse et l’efficacité de 
notre approche, nous présentons plusieurs exemples numériques tout au long de la thése. 
Ces exemples démontrent comment nos solutions numériques résolvent avec précision 
les équations données, ce qui confirme davantage |’efficacité de la méthode que nous 


proposons. 


Mots-clés: Polynémes orthogonaux. Equations intégro-différentielles de Fredholm. 


Méthode de Galerkin. Technique de colocation. Equation non linéaire. 
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Notations 


Vector space. 
Banach space. 
Scalar product. 
Hilbert space. 
Scalar product. 


Norm of z. 


Polynomials of degree less than or equal n. 


Partial differential. 

Lebesgue space. 

Sobolev space. 

Set of complex numbers. 

Set of real numbers. 

Set of integer numbers. 

Set of positive integer numbers. 


Set of strict positive integer numbers. 
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Introduction 


An integro-differential equation is a type of differential equation that involves both derivatives 
and integrals of the unknown function. Such equations are often used in physics, engineering, 
and applied mathematics to model a wide range of phenomena, including heat transfer, fluid 


mechanics, and signal processing, ... (see 6) [36]). 


Solving integro-differential equations can be challenging, as they often do not have ex- 
plicit solutions. Instead, numerical methods are used to approximate the solutions. Among 
these, we can enumerate finite difference methods or spectral methods: Laplace decomposi- 
tion method |5], Legendre-Galerkin method [[8], Bernoulli polynomials [8] , Pseudospectral 
methods, Piecewise linear approximation, Polynomial approximation, Rational approxima- 
tion 2IJ, B-spline method [16] 28], Euler matrix method [25], Exponential spline method 
[22], CAS wavelet [85], Differential transformation [I3], Schauder bases [7], Homotopy per- 
turbation method |89], collocation method [BI], Haar wavelet bases |I7]. 


On the other hand, there are many publications that give the numerical solution by using 
orthogonal polynomials; we find for instance: Chebyshev polynomials [14] [33], Legender 
polynomials [87] 38], Hermite polynomials [24], Laguerre polynomials [80} [32], Jacobi polyno- 
mials [29], new class of orthogonal polynomials [3], Comparison of the orthogonals polynomial 
. However, applying orthogonal polynomials to solve nonlinear integro-differential equa- 
tions can be more complex than solving linear equations. In fact, there are still techniques 
that can be used to approximate the solution using orthogonal polynomials. One approach 
is to use a Galerkin method, which involves approximating the solution as a linear combina- 
tion of a finite number of orthogonal polynomials, and then projecting the original nonlinear 
integro-differential equation onto the space of these polynomials. This yields a system of 
nonlinear algebraic equations for the coefficients of the polynomial expansion, which can 


then be solved using numerical techniques such as Newton’s method or fixed-point iteration. 


CONTENTS 


Another approach is to use a collocation method, which involves evaluating the integro- 
differential equation at a set of collocation points, and approximating the solution as a linear 
combination of orthogonal polynomials that satisfy the equation at these points. This yields 


a system of nonlinear algebraic equations that can also be solved numerically. 


In this thesis, we undertake a numerical research for several forms of nonlinear integro- 
differential equations in which the unknown function and its derivatives are in the nonlinear 
kernel. The current study relies on the use of several kinds of orthogonal polynomials to get 


the required numerical solution. 


The outline of the thesis is as follows: The first chapter aims to introduce some funda- 
mental definitions and theorems in functional analysis, as well as the orthogonal polynomials, 
including their properties that are necessary for the subsequent chapters. 

In the second chapter, we examine regular integro-differential equations with nonlinear 
kernels. In the first section, using Legender polynomials, we give the numerical solution for 


the following equation 


1 
ule) =f) + f K@uue)wW)dy, Ve € [0.1 
0 
u(0) =p 

In the second section, we develop the previous method presented in the first section by 


utilizing Chebyshev polynomials to solve the following equation 


1 


K (z,y, 0(y), vy), v"(y)) dy, Vz € [0,1], 
~(0)=a, (0) =8. 


Through the third chapter, we investigate nonlinear integro-differentials equations with 


s- 
< 
| 
SY 
& 
fe 
o—~ 


weakly singular kernels which takes the following form 


u(z) = f(z) +f p(|z — t|)F (z,t, v(é), o'(t)) dt, z € [a, 4, 
v(a) = U0, 


where 


lim p(x) = +00. 


x20 


10 


CONTENTS 


In the first section, we apply the Galerkin method along with Chebyshev polynomials of 
the second kind to approximate the solution of the nonlinear integro-differential equation. In 
the second section, we use the collocation approach based on Laguerre polynomials to solve 
the equation. 

Finally, in the fourth chapter, we utilize Hermite polynomials to obtain an approximation 


for the solution of the fractional integro-differential equation of the form 


220d —s) + / K (z,t,u(t),D°u(t)) dt, (0) =0, 


where D® denotes the Caputo- Fabrizio derivative of order a [20]. 


During the period of the thesis study, we were able to publish the following articles: 


1. Henka, Y., Lemita, S., Aissaoui, M. Z. (2023). Hermite wavelets collocation method 
for solving a Fredholm integro-differential equation with fractional Caputo-Fabrizio 


derivative. Proyecciones (Antofagasta), 42(4), 917-930. 


2. Henka, Y., Lemita, S., Aissaoui, M. Z. (2022). Numerical study for a second order Fred- 
holm integro-differential equation by applying Galerkin-Chebyshev-wavelets method. 
Journal of Applied Mathematics and Computational Mechanics, 21(4), 28-39. 
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Chapter 1 


Preliminaries 


1.1 General notions 


1.1 General notions 


In order to provide a comprehensive understanding of the material presented in this thesis, 
it is important to establish a solid background in the relevant mathematical concepts. This 
section introduces the reader to some fundamental definitions and theorems that will be used 


throughout the rest of the work. 


Definition 1.1. A norm on a vector space V is a mapping || - || : YV > [0,00) such that for 


allu,v € V and for alla € C: 
(i) \ju+v|| < lull] + lol] (Triangle Inequality), 
(it) \lau|| = |al||u|| (Scalar Property), 


(iti) ||ul]| = 0 if and only if u = 0. 


The pair (V, || - ||) is then called a normed space. 
Definition 1.2. A space B is said to be a Banach space if (B,|| - ||) 1s a complete normed 
space. 


Definition 1.3. Let V be a vector space over C. A scalar (or inner) product is a mapping 
(.,.)fromV x V into C such that 
(i) (u,v) =(v,u) YWuvey, 
(it) (u+w,v) = (u,v) + (w,v) and (Au, v) = A(u,v) Vu,v Ee V,AEC, 
(ii) VueVvV (uu) >0 and (u,u) =0 of and only if u=0. 


We call a linear space with an inner product an inner product space or pre-Hilbert space. 


Proposition 1.1. Any inner product (.,.) defines a norm on V by setting: 


|u|] =f (u,u) Vue V. 


Exemple 1.1. A Hilbert space is a complete inner product space. 
Exemple 1.2. The space H = C” is a Hilbert space with the inner product: 
OD So oti Vu,weC’, 
i=1 
where O: = (015 0i.02.50_) ond = (Uy, toy «sty )- 
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1.1 General notions 


Exemple 1.3. Let H = L?({a,b]). Then H is a Hilbert space equipped with the inner product: 


pe i, f(a)glajdx Yf,g € D((a,8)). 


Exemple 1.4. Let 


H=C/= " = (tn) nens CC, S° Jui|? < ~| : 


i=1 
and the scalar product 
(tu, 0) = So uti Vu,v € £?. 
i=1 
Then H is a Hilbert space. 
Theorem 1.1. (Cauchy-Schwarz inequality) Let H be a Hilbert space. Then 


(v,u)| < |jullflel]| for all u,u € H. 


Definition 1.4. Let M = {u; © H;ie¢ I CN} and (u;,u;) = 0 for all i,j € I with i 4 j. 
Then, M is said to be an orthogonal system. Additionally, if M is orthogonal and ||u;|| = 1 


for alli € I, then M 1s called an orthonormal system. 


Lemma 1.1. (Bessel’s inequality ) Let H be a Hilbert space. If {e; €H,i€ I CN} is an 


orthonormal basis, then for allu © H 


Yd Mexs up|? S lle. 
k>0 
Definition 1.5. A set C ws convex if and only if 
Vuvec, Vee [0,1] du+(1-OVEC. 


Exemple 1.5. In a normed space E, any subspace is a convex set. 


Theorem 1.2. Let K be a non-empty closed convex subset of a Hilbert space H. For any 


element u € H, there exists a unique point v € K such that 
lu — v|| = inf ||u — wl. 
wek 
Theorem 1.3. Let (E,||-||) be a normed space. A mapping T : E > E is called a contraction 
on E if there exists a positive constant p <1 such that 
||T'(u) — T'(v)|| < pllu—v|| for all u,v € E. 
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1.2 Orthogonal polynomials 


Theorem 1.4. [// (Banach’s Fixed Point Theorem). Let (E,|| - ||) be a Banach space 
and let T : E — E be a contraction on E. Then T has a unique fixed point x € E, 1.e. 
1) =a 


Definition 1.6. [70/ (Caputo-Fabrizio derivative) Let a be a real number from the open 
interval (0,1). The fractional Caputo-Fabrizio derivative of order a for a function u belonging 


to the space H'(0, 1] is as follows: 


Du(z) = 2 : - is axe | SW (2- 5) wide, 


Definition 1.7. 26/ (Caputo-Fabrizio integral) Let a be a real number from the open 


interval (0,1). We define the fractional integral of order a using the Caputo-Fabrizio operator 


for a function u belonging to the Sobolev space H'(0,1] as follows: 


T°u(z) = (1-a)u(z) +a als). (1.1) 


Lemma 1.2. Given a real number a such that 0 < a < 1 and a function u in the Sobolev 


space H'*(|a,b]), the following identities hold: 


T* (D®u(z)) = uz) — ula), 


D® (L“u(z)) = u(z) — exp c i . 


Proof. See 


1.2 Orthogonal polynomials 


In this section, we address the topic of orthogonal polynomials, which are a special type of 
polynomials with the property that they are orthogonal with respect to a particular inner 
product. We explore some of the basic properties of orthogonal polynomials, including their 


recurrence relations and differential equations. 
Definition 1.8. The two functions u and v are said to be orthogonal with respect to the 
weight function w(z) on |a, b] 
b 
MUO) = / w(z)u(z)olz)dz = 0. 


Proposition 1.2. Let (Pi)iso be a family of orthogonal polynomials. Then, there is a recur- 


rence relation between Prii, Py, and Py-1 


Yn €N*, Fan, bn, Cn €R, Pasi = (GnX + bn) Pa + CnPr-1- (1.2) 
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1.2 Orthogonal polynomials 


Proof. The family (X P,,, Pr, Pn1,-.-, Po) is a family of polynomials having different degrees 


and so, it is a free family of R,41[X] with n + 2 vectors. We deduce that this family is a 


basis of R,4,[X] and there exist reals ap, bn, Cp and a; for 0 <i <n -— 2, such that 


n-2 


Paid = Gy X Pp, or ore ae ar Cala a So aiP,. 


i=1 


We use the orthogonality of (Pi)iso 
VWO<i<n-2, <Piui,P,>=a, < XP,,P, > +a;||P |? =0. 


According to the expression of the scalar product < XP,,P;) = (P,,XP;), since XP; € 
Rn-1[X] we have (XP, P;) = 0, then 


YO<i<n-2, a, =0. 


1.2.1 Legendre polynomials 


The Legendre polynomials take the following form by the Rodrigués formula: 
iS 9 n 
Theorem 1.5. The Legendre polynomials L,(z) are given by the following recurrence formula 


(n+ 1) Ln4i(z) = (2n4+ 1)zLn(z) —nDyi(z) In(z)=2, Lo(z) =1. 


Exemple 1.6. With n = 4 we get the following six polynomials: 


Lo(z) =1, 

Iy(z) =z, 

Lo(z) = ; (32? — 1), 

L3(z) = : (52° — 32), 

Lge) : (352* — 302? + 3), 
L3(z) = : (632° — 702° + 15z). 
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1.2 Orthogonal polynomials 


Legendre Polynomials Legendre Polynomials 


T,00 
° 
T(x) 


ie 05 Se 


4 -0.5 0 0.5 4 1 05 0 0.5 1 


x x 
Legendre Polynomials Legendre Polynomials 
f 


x) 


Figure 1.1: Legendre polynomials for 1 <n < 4. 


Proposition 1.3. Legendre polynomials are orthogonal on |—1,1] and 


: ay =), 
fl Lz) Ej(z)\dz =< 2+1 
-1 0, ifs. 


Proof. For i<j 


fs len Slew)e [Eile Biles 


1 
Since 
d' 9 i) @! 2 j : . dk 2 j 
Ta (2 —1) al (2 —1) ] = 0 because Vk € (0,1,...,3-1), a5 (2 —1) 


then the relation(1.4)) becomes 


[ 2 (2? _ 1)'| & (2? = 1)'| dz = -f st (2? = 1)'| ee (2? 7 1)"] Ze. 


1.2 Orthogonal polynomials 


By integrating in the right-hand member 7 times, we get 
l a di j fg aod j 
Lele gle-v}e-cy [ale ale] 
qits-1 4 de ait 
+fmale-9] Ble-]], 


-(-1) [ [2-1] (=a 


1 


Since ae 
_ (2? - 1)'] Qi j, 
then, 
[ele] gle-o]e-0 
Ifi=j, 


(n+ 1) I. En4i(z)Ln-1(z)dz — (2n + 1) ie 2Dn(z)Ln—1(z)dz + nf L?_,(z)dz = 0. (1.6) 
(n+1) I. Engi (2) Ln(z)dz — (2n + 1) I, zL?(z)dz+n I. Lyi (2) Ly 2) dz = 0. 
2Ly(z)Lnyi(z)dz+n I. Lys 2) Lina (2)d2 =O. (127) 


1 


nt) f thalerde— n+) | 


—1 
Since 
1 1 
/ Enyi(z)En-1(z) dz = 0, f Ling1(2) Ln (z)dz = 0, 
4: -1 


1 


[ Ln-1(z) En (z)dz = 0 and i To ig de = 0, 


1 1 
then from (1.6) we deduce 


1 


—(2n + 1) is 2L,(2)En_-1(z)dz + nf i? _,(2)dz=0, 


1 


Therefore, 


and from (1.7), we deduce 


1 


nf L? (z)dz — (2n — 1) / Z1m—1(2)Ln(z)dz = 0. 


1 -1 


Hence, 


1.2 Orthogonal polynomials 


Thus 


By recurrence, it comes 


1 7 l 1 7 9 
I. i (z)az= (Qn +1) +) i. Ly lz)\dz = Qn+1) +1) 


Corollary 1.1. The Legendre polynomials can be defined by a differential equation for any 


integer n, such that Ly(z) is a solution for the differential equation: 
(1 — 2°) P” —2zP’+n(n+1)P=0. 


Proof. We have 
Li, (z) = 2nz (2? - iy : 


after multiplication by (2? — 1), we obtain 
(2? — 1) Li (z) = 2nzL,,(z). (1.8) 


By differentiating the equality (1.8) n+ 1 times, with Leibniz’s formula, we find 


(2 — 1) 2 foale)) + (n+ N22) Se thal] + Ye) iz, (2) 
qrtl qd” 
=2n | A [nll] + (01) Leal) 
(| F eatal) +22[F esl] -ntn+ nF tent =o, 


1 
that gives, after multiplication by a ion 
ni2” 


Vn EN, (1 — 27) LP (z) — 22D}, (z) + n(n + 1IL,(z) = 0. 


1.2.2 The first kind of Chebyshev polynomials 


The Chebyshev polynomials of the first kind are given by the Rodrigués representation 


—1)"Vn(1— 22)? an n-1/2 
Utz) at Lived ee (1 - 2?) lp (1.9) 
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1.2 Orthogonal polynomials 


Chebyshev Polynomials of the First Kind 


T(x) 
g 


4 0.5 0 0.5 1 


x 
Chebyshev Polynomials of the First Kind 


T3() 
° 


Chebyshev Polynomials of the First Kind 


cal 0.5 0 0.5 1 


x 
Chebyshev Polynomials of the First Kind 


Figure 1.2: Chebyshev polynomials for 1 <n < 4. 


Definition 1.9. The first kind of Chebyshev polynomials of order n is given by 


U,,(z) = cos (ncos~*(z)) , 


Vz € [-1, 1]. 


Theorem 1.6. The Chebyshev polynomilas of first kind can be defined by the following re- 


currence relation 


Unsi(z) = 22U,(z) — Un-i(2), 


Uo(z) = 1, 


for n2 1, 


Proposition 1.4. These polynomials are orthogonal with respect to the weight function 


w(z) =1/V1— 2? on the interval [—1, 1], and 


Proof. We have 


w/2, nin F 0, 


0, eM. 


cos(né + m@) + cos(n@ — m@) = 2cos(n@) cos(mé). 


1.2 Orthogonal polynomials 


We put z = cos@. Then 


dz = —sin@d@ 
dz 
—dz = sin 0d = \/1 — (cos@)?d0 => dd = — Jos 
a Tn(z)Tm(Z) -[% Bee - cos(n@) cos(m0)dé 
3 === ‘fie 0 
Then ifn Am 
| cos(n@) cos(m@)dé = 0, 
0 
LHS 20 


| cos(n@)?dé = i wel conn) 9 = & 


i cos(n6)*d0 = a) 
0 0 


Corollary 1.2. The Chebyshev polynomial of degree n is a solution for the differential equa- 


ifn=m=0 


tion: 


(1 — 2?) P"(z) — zP"(z)+-n?P(z) =0, Wz e [-1,1]. 


Proof. Let n € N. By deriving the equality T;,(cos 0) = cos(n@), we get 


VYOER, (—sin6)T’(cos@) = —nsin(n6). 


By deriving this equality a second time 


V8ER, (—cosé)T"(cos 6) + sin? OT” (cos 0) = —n? cos(nb) = —n?T,,(cos 6), 
Vz€[-1,Y, (-2) T(z) + 1-2) Biz) = —n7T,,(2z). 


Hence, 


Wn EN, Vz€ (0,1, (1-27) T(z) —2Ti(z) +n7T,(z) = 0. 


1.2.3. Chebyshev polynomials of the second kind 


The Chebyshev polynomials of the second kind U,,(z) are given by the following formula: 


(—1)"(n + 1) Vr d” (1 _ zene 


U,(z) = 
(z) gn+1 (n zh i)! (1 = 22)? dzm 
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1.2 Orthogonal polynomials 


Theorem 1.7. The second kind of Chebyshev polynomials are defined by the following recur- 


rence relation: 


Un+i(z) = 22U,,(z) — Un-1(2), 
Uo(z) = 1, U;(z) = 2z2,n > 1. 


(1.10) 


Exemple 1.7. The first siz Chebyshev polynomials of the second kind are 


Uo(z) = 1, 

Us) = 22 

U,(z) = 427 — 1, 

Us3(z) = 82° — 4z 

U,(z) = 16z* — 1227 +1 
Us(z) = 322° — 322° + 6z. 


Proposition 1.5. The present polynomials are orthogonal with respect to the weight function 


w(z) = V1 — 2? on the interval [—1, 1]. Moreover, we have 

Tse 

. 5 LF): 
/ w(t)U;(t)U;(t)dt = 4 2 

= 0, iFf. 


Corollary 1.3. The Chebyshev polynomials of the second kind of the degree n is a solution 


for the differential equation: 
(1 — 27) y” — 3zy' + n(n + 2)y =0. (1.11) 
Proof. By setting z = cos(@), we have 
U,,(cos(@)) sin(@) = sin((n + 1)), 
and by differentiating this equality twice, we get 
—sin(@)U,,(cos(9)) — 3 cos(@) sin(@)U"(cos(@)) + sin(@)°U/"(cos(8)) = —(n + 1)? sin((n + 1)8). 
Thus, 
—3.cos(@) sin(@)U"(cos(@)) + sin(@)?U“"(cos(@)) = —n(n + 2) sin((n + 1)8), 
and by dividing by sin(@), we find 
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2 Chebyshev Polynomials of second kind 3 Chebyshev Polynomials of second kind 
\ / 
15 25} \ / 
1 2t\ 
05 15 
oe = 
~~ oOo i | 
| a 
0.5 0.5 
4 0 
15 -05 
2'- a = 
1 -0.5 0 0.5 1 4 0.5 0 0.5 1 
x x 
a Chebyshev Polynomials of second kind % Chebyshev Polynomials of second kind 
f | 
| 
3 ] all | 
/ \ | 
2 
3 | | 
\ 
1 ieee \ 
2} | 
= 0 — \ i 
=) 
1 \ = 
i! — \ 
0 
2b] \ 
/ 
3t/ Af 7 
/ 
4 2 
-1 0.5 0 0.5 1 4 0.5 0 0.5 a 


V8 ER —3cos(8)U’(cos(6)) + sin(0)?U"(cos(0)) = —n(n 4+ 2) 


Then, 


Vz €[-1,l] (1-27) UN (2) — 32U}(z) + n(n + 2)U,(z) = 0. 


1.2.4 Hermite polynomials 


The Hermite polynomials can be defined by the Rodrigués formula 


Theorem 1.8. The Hermite polynomials H,,(z) satisfy the following relation: 


22Hn+1(2) — An+2(Z) “Ir 2(n at 1)H,(z). 
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(1.12) 


(1.13) 


1.2 Orthogonal polynomials 


Hermite Polynomials 


4 0.5 0 0.5 1 


x 
Hermite Polynomials 


Hermite Polynomials 


0.5 


0 0.5 1 
X 
Hermite Polynomials 


Figure 1.4: Hermite polynomials for 1 <n < 4. 


Exemple 1.8. The first six Hermite polynomials H,,(z) are 


? 

zy 
—? 
e107, 


Proposition 1.6. These polynomials are orthogonal on the interval R with respect to the 


weight function w(z) =e * 


Lemma 1.3. We have 


(oe) 


Proof. The exponential generating function is given by 


2 
. Moreover we have 


-++0o 
/ exp (—u’) du 
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1.2 Orthogonal polynomials 


. An (z)w” 
G(z, w) = exp (2zw — w’) = d. an, 
We have 
+0o 
/ G(z, w)G(z, t) exp (—2z”) dz = V7 exp(2ut), 
and then, 
+oo ¢m w ; 
S- ee) S° sy inl?) exp (—27) dz = V7 exp(2ut). 
~% m>0 n>0 
Thus, 


2 (= — [- H(z) Hm(z) exp (—z”) ts) w” = Vn exp(2wt)dz 
=Vr>> (Fr) w", 


n>=0 


and therefore, 
» — H,(2)Hm(z) exp (—27) dz = x24". 
Hence, 
+00 


H,(2)Hm(z) exp (—27) dz = 0, ifn 4 m. 


—co 


On the other hand, if n = m by multiplying the equality (1.13) respectively by H,,(z) and 
H,-1(z) we deduce 


H?(z) — 2zHpn_1(z)Hn(z) + 2(n — 1) Hn-2(z) H(z) = 0,n > 2, 


and 
Asi (2)Ap-a(2) = 22H, (2)8, (2) + 2n BZ, (¢) =0,0> 2 
Therefore, 
H?(z) + 2(n — 1)Hn_o(z)Hn(2z) = Hnoi(z)Hp_1(z) + 2nH?_,(z). (1.14) 


Multiplying by exp (—z7), we find 
H?(z) exp (—27) + 2(n — 1) exp (—z”) An—2(z) An (z) 
— exp (—27) Hnyi(z)Hn-1(z) — 2n exp (—z”) H?_,(z) 
=0. 
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1.2 Orthogonal polynomials 


The integration of the latter gives 


+00 


/ * exp (—27) Ha(z)dz + 2(n — 1) / exp (—27) Hn—2(z)Hn(z)dz 


(oe) —co 


+00 


— / ~ exp (—z’) An41(z)An-1(z)dz — 2n | exp (—27) He i (z)dez 


oe) —co 


=0. 


Due to orthogonality, we have 


+00 


/ * exp (—2*) Hn2(z)Hn(z)dz = 0 et / exp (—2*) Hn4i(2)Hn—1(z)dz = 0. 


Therefore, 
+00 +o0 
/ exp (—z”) H?(z)dz = 2n | exp (—2*) H2_1(z)dz. 


By applying this formula n — 1 times, we get 


me H?(z) exp (—z7) = 2r-1n! [- exp (—z’) H?(z)dz =2""n! f exp (—z’) (2z)?dz. 
= > ~ (1.15) 
But, 
+00 oo 1 “00 ‘too 
afe exp (—2”) (z)*dz = afe zexp (—27) adz =4 5: exp (-2)] : + af exp (—z*) dz 


By substituting this result in (1.15) , we obtain 


0° 
H2(z) exp (—27) dz = V/n2"nl. 


—co 


Corollary 1.4. The Hermite polynomial of degree nH,,(z) is a solution for the differential 
equation: 


y” — 2zy' + 2ny = 0. 


Proof. We have 


d” _.» 
H,, = f(=1)\* 2 —z 
(2) = (-1Y"e" Ze 
and then, ; 
/ n 2 d” —2? nm 2 d’t —2? 
H(z) = (—1)"2ze aa +(-1)"e qe 
Multiplying by a gives 
/ —2? n d” —2? n oe —2z? 


1.2 Orthogonal polynomials 


Differentiating (1.16) leads to 


dq” qvti 
22H! (z)e"* + H"(z)e"* = (1) a 
VAL 


and using the Rodrigués formula (1.12), we find 


—22H) (z) + Ai(z) = 2H, (2) — 22 Angi + Anse. 


Therefore, by the recurrence relation (1.13), 


H(z) — 2zH,,(z) + 2nH,,(z) = 0. 


1.2.5 Laguerre polynomials 


The Rodrigués representation for Laguerre polynomials is 


Theorem 1.9. The Laguerre polynomials satisfy the recurrence relation: 
(n+ 2)Dnyo(z) = (2n+ 3 — z) Lagi (z) —n2+1L,(z). 


Exemple 1.9. The first six Laguerre polynomials are 


1 
La@y= - (—2° + 92? — 1824+ 6), 
1 


La(2) = 5y (2° — 162" + 722” — 962 + 24), 
1 
Ls(2) = 59 (—z° + 252* — 200z* + 6002? — 600z + 120) . 


Proposition 1.7. The Laguerre polynomials are orthogonal with repect to the weight function 


w(z) = exp(—z) on the interval [0, +00. Moreover, we have 


1.2 Orthogonal polynomials 


Laguerre Polynomials Laguerre Polynomials 


4 35 
27 30 
25 
oF 
20 
2 
= ioe 
= se 18 
= al 
4 
10 
6} 
5 
8 0 
10 5 
2 0 2 4 6 8 10 2 1) 2 4 6 8 10 
x x 
Laguerre Polynomials Laguerre Polynomials 


L() 
ul 
= 


Figure 1.5: Laguerre polynomials for 1 <n < 4. 
Corollary 1.5. The Laguerre polynomial L,,(z) satisfies the following differential equation: 
zy’ + (1—z)y’+ ny =0. 


Proof. The generating function is 


g(2,t) = exp i | = So Lalz)2", 


n>=0 


and we have 


O O 
(1-1 2 (z,t) = (1-2—dglz,8), (€- 1 ZE(e,4) = tel, 9), 
Og Og 
772 ) ia oe) “7 
Then, ; 
O"g Og Og _ 
aaa (2,t) + (1 a) 5 (2, t) + ta (2, t) = 
Therefore, 
S> (2ER(z) + (1 — 2) Lh, (2))t? +45 nLa(2)t” = 0 
n>0 n>1 
Hence, 


28 


1.2 Orthogonal polynomials 


29 


Chapter 2 


Regular Nonlinear Fredholm Integro 


Differential Equations 


30 


2.1 Fredholm integro-differential equations of the first order 


2.1 Fredholm integro-differential equations of the first or- 
der 


In this section, we present both analytical and numerical investigations of Fredholm integro- 


differential equations of this particular type: 


ule) =f) + f Keuul)w) de, (2.1) 
0) =p, 


u( 


where K, “ € C((0, 1]? x R*), f(z) € H"((0,1]) and u(z) € H'[0, 1]. This section aims 
to investigate the existence and uniqueness of the solution of problem (2.1)) and propose 
a numerical process for its approximation. To achieve this goal, we analyze the necessary 
conditions for the existence and uniqueness of the solution. Then, we apply our numerical 


method that can approximate the solution with a high level of accuracy. 
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2.1 Fredholm integro-differential equations of the first order 


2.1.1 Existence and uniqueness 


Consider the Sobolev space H = H'(|a, bj, R) equipped with the following norm 


Vu € [feelin = Mell z2tay) + Mell c2j0.9)- 


As a starting point, the following assumptions are made: 


a,, 6, > 0, such that r = 0,1Vz,x € [0,1], for all u,v, u,v € R, 


|K(z, 2, u,v) ~~ K(z,2, 4%, d)| < au |u _ a + Bolv ~— vl, 


|O.K(z,2,u,v) — 0,K(z,2,t,d)| < aylu — a] + B,|v — a], 


0<y=max {aot a1, 8+ 6} <1. 


Let f € H and define the operator: 


¥z€ [0,1], L:w#—H 


wr Leyla) = 42) + f KGy.utu) alo) dy 


with [B(w)(2) = F@) +f a.k (2.9. uly), 4") dy. 


Theorem 2.1. The equation (2.1) admits a unique solution in H based on the assumption 
a: 


Proof. Let vy, € H. Then, for all z,x € [0, 1] 
|K (z, 2, p(2), p(x) — K (z, 2, 0(2), Y'(x))| < aolp(x) — ¥(x)| + Boly'(x) — ¥'(@)|. 


It follows that for all z € [0,1] and by Cauchy-Schwarz inequality 


llL(y)(z) — LW)(2)| < aole(y) — YM) + Foley) — YO! 
€ aglly — Vilz2 + Bolle’ — W'llz2 - 
IL(y)(z) — L@)(z)? < aglle — WIliz + 66 lle! — Ullz2 + 2ooBolle — Viz lle’ — ¥'llze 
Le) — LW)Ilz2 < (aolle - Vllzz + Bolle’ - Ulla)” - 
Wenee: 
L(y) — LW) Ilzz < aolle — Pilz + Bolly’ — Y'Iiz2 - 


Similarly, we can find 


IL(¢)' — Lb) 'IIz2 < aalle — vile + Ally - ¥'llze- 
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2.1 Fredholm integro-differential equations of the first order 


and then, 
|EZ(~) — L@) ln < ile — Vile. 


Since y < 1, the equation (2.1) has a unique solution. 


2.1.2 Numerical study 


We obtain the numerical solution using Legendre wavelets and the operational matrix of 
integration. The Galerkin method is applied to obtain a nonlinear algebraic system, which is 
subsequently solved using the iterative method. Additionally, we provide several illustrative 
examples. 


The Legendre wavelets can be defined as: 


ba i i 
9's ./25 FIL, az —2i+1)), ag Te 
big(2) = i ))> pi 2k (2.2) 
0, otherwise 


where k € N*,j € NO < j < m—1,i = 1,2,...,2°7', and L;(z) is the Legendre 
polynomial with degree j. Furthermore, this family of Legendre wavelets {¢;,;} also defines 


an orthonormal basis for L?({0, 1)). 


Function approximation 


Using the Legendre wavelets basis, each function u(z) in L?((0,1]) can be expressed by the 


following formula: 


where cj; = (u, £;,;), and (.,.) denotes the inner product in L7((0, 1]). 
By truncating the infinite series (2.3), we can write: 
2k 1 m—1 
Um 2) ~~ a S- Cigligl2) = CT Pz), (2.4) 


i=1 j=0 


where C7 and P,,(z) are 2*~'m x 1 matrices given by: 


T 
GC = [e10, C115 +++ 5 C1ym—15 €1,0) C11) +++) Clm-1)--- > Cok-1 9, 2 ieee tea 5 


T 
Pr (z) => [é1.0, C14, sae i fignals Lids C14, sae si ews sae ; for=i.0, Cok—1,1 sae , Cok_1.m—1| 7 
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Figure 2.1: Legendre wavelets for k=1 and 1<n<4 


Theorem 2.2. [{0/ Let u € L7[0,1] be with a bounded second derivative, (e.g., u"(z) < M). 
The truncated series Um(z) given in (2.4) converges uniformly to the function u(z). Moreover, 
242M 
cig] S I 
2°#+1(25 + 1)(27 — 1) (29 — 8)2 


Operational matrix of integration 

As a particular case, we take k = 1 in (2.4), to get: 
ct = [co, Cly eens Cini | Fi 

Pa) = lol 2), G(2) 2005 La A) | 


The matrix M/,, is defined as a matrix consisting of the coefficients of Legendre wavelets. 
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2.1 Fredholm integro-differential equations of the first order 


Pye Ve: os (ane gaa 
0 2V3 -6V5 
0 0 65 
Pe NA — 0 
Oi <0 - uel oe 
((m — 1)!)? 
In the canonical polynomial basis X,,(z) = (1, a eee aay the vector P,,(z) of Leg- 


endre wavelets can be rewritten as follows: 


Consider the matrix N which contains the coefficients of the integral of canonical poly- 


nomial basis: 


O10 0. + 0 
0 0 : 0 0 
a 1 
N=]000 = - 0 
3 
1 
00 0 0 — 
m 
Then, we have the following integration matrix : 
i. CT P,(s) = CT N Pryi(z) = C7 Qm(z), (2.5) 
0 


where Qm(z) = NPn4i(z) and N is given by: 


N = M3 NMn4i. 


Method description 


Let be given the integro-differential equation: 


ul) =f) + f Kevu)wo) dn WO) =e 


First, let us differentiate the above equation with respect to the variable z to get: 


wle)= f+ f a.K @..uly).wW) dy (2.6) 
3 
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2.1 Fredholm integro-differential equations of the first order 


Now, we approximate the unknown function u'(z) of equation (2.6) by using Legendre 


wavelets (2.4) : 
i bi OC Pot): (2.7) 


Next, we integrate the equation (2.7) with respect to the variable ¢ from 0 to z : 
Um(2) = @+ C*Qm(z), (2.8) 


where Q,,(z) is described in (2.5). By substituting relations (2.7) and (2.8) into the 
equation (2.6)) we obtain the following non linear equation corresponding to the m unknown 


coefficients (c;,2 = 0,1,...,m—1): 


1 
CT Pn(z) = f'(2) +f OK (z,y,0+ C* Qm(y), C7 Pm(y)) dy. 
0 
To obtain the coefficients C? = [Co, C1,---;m—i], we apply the Galerkin projection method 
corresponding to the Legendre’s wavelet basis, i.e. by multiplying the above equation by 
é;(z),j =0,1,...,m— 1 and integrating with respect to the variable z from 0 to 1 , we get 


the following non linear algebraic system: 


Co -| f'(z)bo(z)dz +f / £o(z)O.K (z,y, 0+ C7 Qm(y), C7 Pm(y)) dydz, 


“1 = / f'(2)la(z)de + / / HK (not CTQa CPP dude, 


1 1 1 
oe / f'(2)lm-a(2)dz + i / bd OK Gio COx) CP) aude. 


The exact solution C7 can be difficult or even impossible to obtain in several cases. 
Therefore, approximative methods are widely used, generally iterative methods. In our study, 
we choose the Picard successive approximations method. We start with the initial vector Cy 


and consider the sequence of vectors (Cf) We then consider the following system: 


keEN* 


1 1 1 
cht =f p(eptele)de+ ff bo(2)0.K (2.4.0 + CEQm(u),CPPm(y)) dude 
0 0 0 


1 1 1 
Gf = | f'(z)er(z)dz + : | £1(z)O.K (z,y, 0+ CE Qm(y), Ch Pm(y)) dydz, 
(0) 0 0 


1 1 1 
cht, | Pb iOiles | | ln1(2)02K (z,y, 0+ CLQm(y), CT Ply) dyde. 
: — (2.10) 
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2.1 Fredholm integro-differential equations of the first order 


One can solve the algebraic system (2.10)) by using the Picard successive method to get the 


coefficients (c;,i = 0,1,...,m 


represents an approximate solution u* (z) to the main equation (2.3). 


Convergence analysis 


—1). Then we substitute them into the formula (2.8), which 


Theorem 2.3. The numerical solution Um converges to the true solution u in the Sobolev 


space H. 


Proof. We will discuss the demonstration of this theorem in a similar way through the next 


section. 


In this section we only demonstrate the convergence of the system (2.10). 


Consider the Banach space R™ equipped the norm 


m—-1 
de lel’. 
i=0 


Vac € R™ 


Let’s suppose that: 


Vz,x € [0,1], for all a,a,b,b € 
|OnK (2; x, a,6) — 0,K(z, 2, G, b) 
d=a+P<1. 


Iltll2 = 


R, 


< ala — a| + Bb — 8, 


Theorem 2.4. Under the assumption (A), the nonlinear system (2.9) has a unique solution. 


Proof. Let us define the following operator: 
, 73): R™ —> R” 


oan (0) = f ae 


Let C7? and G” be two vectors from R™. 


T: (T,,To,-+- 


wrod ff 


We have 


IZ (C7) -7(@)I, < 2 CTP)’ - 2 (eT PH)’ 
< a||MC7 P(t) — MG? P(2)||_,. + 6 ||C7 PC) 
< al|M|| ||C7 P(t) — G7 P()||,. + B||C7 PO) 
< aljM|||C” - G7], + ||” - "|, 
<(e+6) |G" = 
<6||c"— "I, 
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2)0.K ( (zy, 


MC P(y),C*P(y)). 


Bessel’s inequality 
—G*P(t)||,2 


— GP (OI 


2.1 Fredholm integro-differential equations of the first order 


where the operator L is defined above. Since 6 < 1, by Banach’s theorem [I-4] | the system 
2.9) has a unique solution. 


Theorem 2.5. Under the assumption A, for any initial vector Cis the sequence (CY) oe 


converges to the vector C’. 


Proof. We have 
[Cena — "|| s 4||Ce — OF]. 


Therefore, by recurrence on k, we get 


[Cees — OF SP |[Cg — C*. 


id — Cc" || — 0 when k — +00. 


Corollary 2.1. |r, — ull, < |e — tli + ||Um — ullz, + 0 when k,m — +00. 


2.1.3. Numerical examples 


Here, several numerical examples are given to demonstrate the efficiency of our proposed 
method. We mention that the numerical results are computed by using the following error 


function: 


2 
Em = ||Um(z) — alley ({0,1]) = 4 llum(2) — 2) INz2 ({o,1]) + Ilia (Z) — u'(z)\lz2((0,1))> 


where u(z) is the exact solution and u,,(z) the approximate solution given by our proposed 
method. The results in the tables below are obtained when \|Coa = Cr || < e, for different 
values of €. 


First example 


Let us consider the Fredholm integro-differential equation: 


f(z) = ze* - - [sin?(1 + z) — sin?(z)]. 


The exact solution is u(z) = ze’. 
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in b2=10" Peal £610" 2,,6= 107 
2 24042e—02 2.4042e—02 2.4042e—02 2.4042e — 02 
3 1.8363e—03 =1.8363e —03 —-1.8363e — 03 1.8363 — 03 
4 1.0855e—04 1.0855e—04 1.0855e—04 1.0855e — 04 
5 5.2218e—06 5.2218e—06 5.2218e—06 5.2457e — 06 
6 2.1146e—07 2.1146e—07 2.1193e—07 5.4274e — 07 
7 7.389le—09 7.3891le—09 1.5907e—08 4.999le — 07 
8 2.3076e—10 2.3117e—10 1.4089e—08 4.9985e — 07 
9 4.4605e—11 4.6703e—11 1.4087e—08 4.9985e — 07 


Table 2.1: Numerical results for the first example . 


Second example 


Consider the following integro-differential equation: 


u(z) = sleos(z +1) +3cos(z)] + i. mi 


«(0) = 


0. 


The exact solution is u(z) = sin(z). 


Third example 


Consider the following equation: 


where 


and the exact solution is u(z) = | sin(27z—7)|. In this example, to obtain an approximate 


solution um(z), we must take k = 2 in the Legendre wavelets formula (2.2), because u’(z) 


sin(z + y) 


1+ u(y)? + u(y 


u(z) = f(z) - / inte) 3 + u?(y) + 0.257-? (w'(y))*dy, x € (0, 1], 


i)\= 5 sin(2) +|sin(27z —7)|, 


has a discontinuous point at z = 0.5. 
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2 ty: z € (0, 1], 


2.1 Fredholm integro-differential equations of the first order 


i. Hips 107°? J e100! ee ea Or® 
2 1.5073e—03 1.5073e—03 1.5073e—03  1.5073e — 03 
3 1.7277e—-04 1.7277e—04 1.7277e—04 1.7277e — 04 
4 4.7584e—06 4.7584e—06 4.7584e—06 4.7588e — 06 
5 3.6343e—07 3.6343e—07 3.6344e—07 3.6913e — 07 
6 7.1237e—09 7.1237e—09 7.3462e—09 6.4994e — 08 
7 4.0789e—10 4.0789e—10 1.8400e—09 6.4604e — 08 
8 8.4572e—12 8.5713e—12 1.7943e—09 6.4602e — 08 
9 5.8232e—12 5.9877e—12 1.7943e—09 6.4602e — 08 


Table 2.2: Numerical results for the second example . 


Ge Bee 10. Bee 0 ee 
3 1.7436e —02 1.7436e—02 1.7436e — 02 
4 1.7436e—02 1.7436e—02 1.7436e — 02 
5 2.6345e—04 2.6345e—04 2.6347e — 04 
6 2.6345e—04 2.6345e—04 2.6347e — 04 
7 2.2965e—06 2.2973e—06  3.9486e — 06 
8 2.2965e—06 2.2973e—06 3.9486e — 06 
9 2.8432e—07 2.9053e—07 3.2246e — 06 


Table 2.3: Numerical results for the third example . 
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Figure 2.3: u’(y) vs ul,,(y) with m = 9 for the first example. 
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Figure 2.4: u(y) VS Um(y) with m = 9 for the second example . 
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Figure 2.5: u(y) vs ul,,(y) with m = 9 for the second example . 
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Numerical Solution 


Exact solution 
* Numerical solution 


0 0.2 0.4 0.6 0.8 1 


Figure 2.6: u(y) vs Um(y) with m = 7 for the third example . 


Exact solution 
* Numerical solution 


Numerical Solution 


Figure 2.7: u/(y) vs u),(y) with m = 7 for the third example . 
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2.1 Fredholm integro-differential equations of the first order 


Interpretation of results: 


Tables 2.1] 2.2] and 2.3} show the error function E,, for different values of m for the previous 
examples. In all cases, it seems that accuracy increases as m increased. On the other hand, 
Figures [4.2}- . , 2.7]show the comparison between the exact solution u(s) (with its derivative 
u’(s) ) and the approximate solution u,,(s) (with its derivative u/,,(s) ) for the three examples. 
It also appears that the exact and approximate solutions are almost identical. So, we confirm 


from these results the efficiency and validity of our proposed method. 
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2.2 Second order Fredholm integro-differential equations 


2.2 Second order Fredholm integro-differential equations 


Our aim through this section is to introduce a numerical study for the following type of 


integro-differential equation of a second order: 


of K (z,y,W(y),¥"(y), 0") dy, 
~(0) =a, = B, 
where K,0.K,02K € C° ((0, = x ,a, 8 €R,v(z), f(z) € H?((0, 1)). 


We utilize Chebyshev wavelets to obtain an approximate solution, and to improve accu- 


(2.11) 


racy, we extend our analysis to wavelets defined on subintervals. We also construct a double 
operational matrix of integration over different subintervals. Finally, we provide several ex- 


amples to demonstrate the efficiency of our approach. 


2.2.1 Chebyshev wavelets 


The Chebyshev wavelets are defined as follows: 


Bie deg i-1 i 
yal BOG), Eheecgh, 
0, otherwise 
where 
1 
p=, 


fork € N*,j7 EN,O <9 <n—-1,i = 2?,p=0,---,k—1. 0; is the Chebyshev polynomial 
of degree j. The Chebyshev wavelets denoted by 6;,; form an orthonormal basis in the Hilbert 
space L*, ([0,1]) with: 


1 
1 
W,4(Z), a oe Saal 
gk-1 _ 
Wok-1 (2) 5A <2 <1, 
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2.2 Second order Fredholm integro-differential equations 
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Figure 2.8: Chebyshev wavelets fork =1and1<n< 4. 


where w;,(z) = w (2*2 —2i+ 1). The graphical representations of the Chebyshev wavelet 
(2.12) 


charts are shown in Figure 2.8 fork = land1l<n< 4. 
Any function u(z) in L%, ([0,1]) can be written as follows: 


b(z) = >> > i 56:5(2), 


i>1 j>0 


where ¢,; = (2, 6;,;), and (.,.) is the scalar product in L?, ([0, 1]). 


By truncating the infinite series (2.12), we approximate the function w(z) as follows: 
2k-1 n-1 
Vnl2) = D2 > c49i,s(2) = CP(2), (2.13) 
i=1 j=0 


where P(z) and C” are 2°~'n x 1 matrices: 
Ce [e10. Cis +++ 5 C1n—1s €1,0) €C2,15 +++ 5 C2.n—1,-+++5 Cok-19)--- Cate | ) 


and 
T 
P(z) = [91,0, O44, z sO teers 92.0, O23, oe loi 2 , 9or-1 9, Bon—1 4 adie pO ont 94 | 
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Figure 2.9: Chebyshev wavelets for k = 2 andl <n< 4. 
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2.2 Second order Fredholm integro-differential equations 


2.2.2. Operational integration matrix 


We take k = 2. Then, P(z) and C” are simplified as follows: 


Cr = [c10, C115 +++ 5C1n—1, €2,0, C21,--- ome ’ 
PZ) = [01,0(z), Giant). eo Orel 2), 02.0(z), 62.1(z), Pear’ ,92.n-1(Z)] _ 


Let W,, be a matrix of Chebyshev wavelet coefficients: 


Fy Ox 


T 
0 0 4 2 
sen eer Cs ae) 
2 | 0 4Vv2 
F,=— . 
T . . = 
0 O re Ar-l1,/9 
and 
Le) = (gee eese hee eee a : 


Py) = [A1,0(Z), 611(2), ee Oi ft Als 62,0(z), 62,1 (2), neds 92n-1(Z)] : 
So, we can write: 


Pal 2) = Zal ZW 


Consider the matrix N, that represents the integral matrix in the canonical basis: 


G, On 
N= a 
O, Gp 
where ; 
—~ 10 0 0 
1 0 : 0 0 
22 2 
G.=| 0 oo 2 0 |, 
3 
1 a 1 
2” n 
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2.2 Second order Fredholm integro-differential equations 


and 
010 0 - 0 
00240 0 
2 

&=|00 o =. 0 

” 3 
1 
n 


Consequently, we find the initial operational matrix of integration: 


| ; CTP, (€)d& = CP My Pasi(z) = C7Qi(2). 


0 


By the same way we get the double operational matrix of integration: 


[ i CTP, (€)dédx _ CO? Mp Mpys1Prso(z) — C7 Q2(z), 
0 0 


where 


M, = Wo N,Waa- 


2.2.3. Method description 


Consider the following problem: 


(2) = fle) + / K (2,£, (6), WOW") a, 
(0) =a, W(0) =8. 


First, we differentiate the equation in (2.14) twice to get the following equation: 


(2) =f") + / AK (2,6, U(6), WO, v" Oa. 


Next, we approximate the unknown function ~"(z) using Chebyshev wavelets: 
" _ nT 
Yn(z) = C° Paz). 
By integrating the equation(2.16) from 0 to z, we find: 
wi (z) = B+C* Qilz). 
Integrating again (2.17) from 0 to z produces 


ti,(z) = a+ Bz + C* Qo(z). 
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(2.14) 


(2.15) 


(2.16) 


(2.17) 


(2.18) 


2.2 Second order Fredholm integro-differential equations 


Then, by substituting Q.16), 2.17) and into (2.15), we obtain: 
1 
CT Pa(z) = f"(z) +f O2K (z,€,a+ BE + C7 Qa(€),B + C7 Qi(E), C7 Pr(E)) dé. (2.19) 
0 


Let us now multiply the equation (2.19) by 0;,;(z)wi2(z) for 7 = 0,1,---,n—1 and 
i = 1,2, and then we integrate with respect to the variable z from 0 to 1. This gives the 


following nonlinear algebraic system: 


1 1 
Cio = Yio + / i, Bio(z)wia(z)O2K (z,€,a+ BE + C7Q2(E), 8 + C7Q1(E), C7 Py (E)) d&dz 


1 1 
Ci = Yi + / / O41 (z)wir(z)O2K (z,€,a + BE + C7Q2(E), B + C7 Qi(E), C7 P(E) dédz 


cant = ina + f / “Gsn-a(2)wia(2)O2K (2,€,0 + BE + C* Qa(E),8 + C7Q(6), CPPA(E)) dz 
sa (2.20) 
where y;,; = (f”,0;,;). By applying the Picard successive approximations method, we 
obtain the vector solution C’ for the system outlined above. We can subsequently 
substitute this solution into to get the numerical solution for the main equation (2.11). 


2.2.4 Convergence analysis 


In order to establish the convergence analysis for the numerical process outlined previously, 


we begin by considering the Sobolev space H = H?((0,1],IR) equipped with the norm 


Wh eH, |Wlle = Ue llz2(o.a) + IY Mhc2p0,4 + WY Mcp, 


Furthermore, let us consider the following additional assumptions: 


4A,, B,,C, > 0, where r = 0,1,2.Vz,t € [0,1], Vz, %,y,9,5,5 ER 
|K(z,t,£,y,8) —K(z,t,Z,9, 3)| < Ao|w — z| + Boly — 9| + Cols — 3|, 

(S) |O.K(z,t, 2, y,8) — OK (z,t, ZY, 3)| < Alu — a] + Bily — 9] + Cils — 3, 
|O2K (z,t,u,y, 8) — 2K (z,t, &, ¥, 3)| < Aglu — a] + Baly — gl] + Cols — 3], 


2 2 2 
<7 mmf ano B. Lot eal e 
r=0 r=0 r=0 


Theorem 2.6. Under assumptions (S), the numerical solution Wy, converges to the exact 


? 


solution w in the Hilbert space H. 
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2.2 Second order Fredholm integro-differential equations 


Proof. Let T be an operator which is defined as a mapping from H to itself as follows: 
T:H OH 
1 
wore) = [ KE EVOWO.0" Oa +f) 
0 


Consequently, we can express the exact solution w to equation (2.11), along with its first 


and second derivatives ~’ and yw", through the following system of equations: 


lz) =TW)(2), 
o'(z) = Te) (2), 
oy" (z) = T"(p)(2). 
Applying the Galerkin projection method by using Chebyshev wavelets, defined in 2-13), 


we can approximate the preceding system as follows: 


It is clear that 
[Wn(z) — V2) = [En Hn) — TY) = [En (Hn) — Tn) + Tn) -— TW) 
Based on the Cauchy-Schwarz inequality and hypotheses (S), we deduce: 


< Ap / IWn(€) — W(6)| dé + Bo i Wh (€) — Wd + Cy / wt(e) — w"(6)| a€ 


< Ao Ilr zs V2.4 + Bo Iv, ~ V'll20,1 + Co Ion a YI ret0 ; 
(2.21) 


Conversely, the research conducted in [2] assumes the convergence of the sequence S,, = 
Tn (Wn) and establishes the following error convergence rate: 
[Tn (Yn) — T Yn)| < O(n") > 0. (2.22) 
So, from inequalities (2.21) and (2.22), we obtain: 
Pn _ Vl r2,0,1 < Ao I/Yn ~ VI lr2,0,1 + Bo Iv ~ Y'lle200,1 + Co Iv ~ V" l20,1 +0 Ca . 
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2.2 Second order Fredholm integro-differential equations 


Similarly, we find: 
IIe — V'llr20.4 <A; len a Vl 12001 + By IIe a Ul eeoa + Cj Ile i Y' Wn200,1 eo (n) ) 


Ion — O' UIcetoay < Ae lle — Pll zeta) + Be Id — YU ceo) + C2 lve — P'Ulep0,4) + O (nr). 


Therefore, 


O(n") = O(n") + O(n") + O(n?) > 0. 
Furthermore, when 0 < y < 1, we can deduce that: 


O(n“) 
Y 


— 0. 
1- 


Yn — Vllee S 


This indeed confirms the convergence of the approximate solution w, to w in the Hilbert 


space H. 


2.2.5 Numerical examples 
First example 


Let be given the integro-differential problem: 


ORIOE: / K (2,£, U6), WO. W"(O) a, 
(0) = (0) =0 


with f(z) = zsin(z) + = sin(=)(21n(2) — 1), and 
Ku 0, 2532) = -5 sin(w) In c + cos(v)u + sin(v)a — ; sin(v)(z — sin(v)) 


and the exact solution is w(z) = zsin(z). 


Second example 


Let be given the Fredholm integro-differential problem: 


vO)=7, ¥(0)=-F, 
with 
(2z — 1) exp(z) + cos(z) [eos (5 — 2cos(2) + cos ( » #E [ 5 
a —(2z — 1)? exp(z) + cos(z) [as 5) — 2.cos(2) tac | ums ft 1 
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2.2 Second order Fredholm integro-differential equations 


09°, 
Exact solution 
0.8 + * Numerical solution 
0.7 
0.6 
c 
o 
3 05+ 
[e) 
” 
w 
iD 
® 
—E 
5 
Zz 


0 0.2 0.4 0.6 0.8 1 


Figure 2.10: w vs w, for the first example with n = 7. 


n 3 4 5 6 7 
E,, 6.988E — 05 | 1.433E — 06 | 4.934E — 08 | 6.723E — 10 | 5.444E — 10 
CPU time 0.109 0.123 0.154 0.194 0.310 


Table 2.4: Numerical results of the first example. 


and the exact solution is: 


1 
7 (22 lrexp(z), 0<2<5, 
W(z) = 4 1 1 2 
2 
——(2z — 1)" exp(z), 5 <—2= 1, 
n 3 4 5 6 7 
En 7.949E — 04 | 2.759E — 05 | 7.617E — 07 | 1.741E — 08 | 8.503E — 09 
CPU time 0.081 0.108 0.129 0.169 0.372 


Table 2.5: Numerical results of the second example. 
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2.2 Second order Fredholm integro-differential equations 


Numerical Solution 


Numerical Solution 


14, 


Exact solution 


12 * Numerical solution 


0 0.2 0.4 0.6 0.8 1 


Figure 2.11: w’ vs wi, for the first example with n = 7. 


Exact solution 
2 * Numerical solution 
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Figure 2.12: w” vs w" for the first example with n = 7. 
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0.35 
Exact solution 
02+ * Numerical solution 
0.1/5 
0 - 


Numerical Solution 
ro) 
RN 
T 


Figure 2.13: w vs w, for the second example, with n = 7. 


Exact solution 
0 * Numerical solution 


Numerical Solution 


Figure 2.14: w’ vs y', for the second example with n = 7. 
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4 7” 

Exact solution 
ot * Numerical solution 
0 


Numerical Solution 
oy 
T 


0 0.2 0.4 0.6 0.8 1 


Figure 2.15: w” vs wi" for the second example with n = 7. 
Discussion 


Based on both the tables and figures, it is evident that the error function diminishes signifi- 
cantly, particularly when the number n is large. This indicates that our suggested approach 
exhibits greater efficacy when dealing with higher degrees of the polynomial (approximate 


solution) n. 
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Chapter 3 


Weakly Singular Nonlinear Fredholm 


Integro-Differential Equations 


In this chapter, our primary objective is to delve into a numerical approximation method 
tailored for solving a class of nonlinear Fredholm integro-differential equations characterized 


by kernels that possess weak singularities: 


(3.1) 


where 


f €C*((a,)), R), F,0,F € C° ([a, b]? x R?) , 
p(t) € W*"(0,b— a), 
lim p (t) = +00. 
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3.1 Galerkin method 


3.1 Galerkin method 


The Galerkin method and the operational matrix of integration are used to compute the 
approximate solution. Based on the features of the second kind of Chebyshev polynomials, 
we present a technique for removing the kernel singularity during the computing process, and 


then we find the desired numerical solution. 


3.1.1 Function approximation 


Consider the Hilbert space H = L?,[—1,1] that is equipped with the inner product (u,v) = 
1 
i: w(€)u(€)u(€)dE with w(f) = /1 — €?, and the orthonormal basis 


1 
2, 
B= {en y/2ee naotent, 
TT 


where @, is the Chebyshev polynomial of the second kind of degree n described in the 
first chapter. Define the following orthogonal projection: 


Weel] Fan F, 


n (3.2) 
ur > P,(v)(€) = ) (u,xa;) o(€) = C7 PO), 


where 


CP = |e, C15 ** 5 Gals = 0,0) (4=05— 0 (3.3) 
and 


P,(&) = [yo(€), ¥1(€), =  Yn(E) . 


P,, is the space of polynomials of degree less than or equals n. We approximate any function 
u(&) of H as: 

u(€) © (Pav) (€) = C7 Pal€). (3.4) 
Theorem 3.1. Let f(€) € L?[—1, 1], satisfying |f”(€)| < L. Then, the series (8.2) converges 
uniformly to f(€). Furthermore, the coefficients in (8.3) satisfy the following inequality: 


AV27rL 
les | < sa aaa Vn. 2 2. 
(n +1)? 


Proof. We have ; 
n= f fe )eale)w(E)de. (3.5) 
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3.1 Galerkin method 


If we set € = cos@ in (3.5), then we get 


= ee f (cos(@)) sin(n + 1)@ sin 6d 
TO 


ee 
=e / f(cos(@)) [cos n@ — cos(n + 2)6]dd, 


that, after integrating twice by parts, gives 


1 . i 
= al f" (cos @)S;,(0)d0, 


where 
5,(0) = sin 0 (Se= 1)@ — sin(n + ) _ sind (Ses ie mn) . 
n n—l n+1 n+2 n+l wn 
Therefore, 
l i MN 
ln] = Jin / fi (os(6)) $08 
= = / F"(cos(8)) (0) 
L Tv 
< al 15,,(0)| d8 
oe ee 1 (— | =| 
V2n [n\n-1 n+l nt+2\n4+1 6 n+3 
2/2 L 
~ (n? + 2n— 3) 
4V/2rL 
(n+ 1)? 


3.1.2 Operational integration matrix 


The operatinal matrix of integration for n = 4 is obtained as follows: 


[ votoae= | a ; 0 0 0 0] PC) 


I 
0 1 0 0 0 | PC), 


3.1 Galerkin method 


Therefore, 
[ct P dt = CMs), 
= 
where 
Q(z) = [po(2), Pilz), +++ » Pnti(Z)], 
and 
1 
1 5 0 0 0 O 
1 
ae 0 - 0 0 O 
4 4 
1 1 1 
Msy6 = 3 6 0 6 0 0 
: 0 : 0 : 0 
i 0 1 ; 1 
5 10 10 


Similarly, we can find the operational matrix of integration M for an arbitrary value n. 


3.1.3 Method description 
By taking the derivative of the equation (3.1), we get the following equation: 
v(z)=7f') +f sign(z — t)p'(|z — ¢t|)F (z,t, v(t), v'(t)) dt 
+ f pet a.F .t,0(0),u'@) at, 


where 


Now, we approximate the derivative of the unknown function v/(z) by using (3.4). Thus, 


gizy=O P(e). (3.7) 


Integrating (3.7) with respect to z from 0 to z gives 


v(z) = 19 + CTMQ(z). (3.8) 
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3.1 Galerkin method 


Now, by substituting (3.8) and (3.7) into (8.6), we get the following equation for the unknown 
vector C7: 


1 


C* P(z) = f'(z)+ I, sign(z — t)p'(|z — t|)F (z,t, vo + C™MQ(t),C* P(t) dt 


1 (3.9) 
+ / p(|z— t|)0.F (z,t, v0 + C7 MQ(t), C’ P(t) dt. 


i 


By multiplying the equation (3.9) by w(z)y;(z) for i = 0,--- ,n, and then integrating with 


respect to the variable z, we find 
co = (f', Yo) rf fw z) sign(z — t)p'(|z — t|)F (z,t, v0 + C’ MQ(t), C” P(t) dt 
‘fmm s|)O.F (z,t, v0 + C7 MQ(t), C’ P(t)) dt, 
ean rf fw 2) sign(z — t)p'(|z —t)F (z,t, v0 + C7MQ(t), C7 P(t)) dt 


+f fix 2)pilz)p(|2 — |) OF (2, t,v0 + CT MQ(t), CT P(t) dt, 


=F ea) f [wo z)sign(z —t)p'(|z -¢t|)F (z,t, v0 + C7 MQ(t),C’ P(t)) dt 


+f 1 z)pn(z)p(|z — s|)O.F (z,t, vo + C™MQ(t),C* P(t) dt. 
— (3.10) 
Analytically, the first integration for each equation in (8.10) exists because p’ € L'({—1, 1}). 
However, the Matlab software cannot run. Hence, we use the integration by parts and the 


following property for the weight function w(z): 
w(—1) = w(1) = 0. 
Then, the equations are simplified as follows: 
= (f', Yo) +f f we p(\z —t)F (z,t, vo + C7 MQ(t), C7 P(t)) dtdz 


Cc, = (f',¢1) +f f we ‘o(|z —t)F (z,t, Uy + C’ MQ(t), C’ P(t)) dtdz aig) 


Cn 


Cae / . / _(00(2)en(2))' nlle = tDF (2.4.09 + C7 MQ), C7 P(O) aide. 
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3.1 Galerkin method 


The solution of the nonlinear equations system (3.11) can be obtained by using Picard suc- 


cessive approximations, which lead to the following system: 


eS ie) +f f coi ‘(jz — tl) F (z,t, vo + Ce MQ(t), CE P(t) dt, 
a = 1,03) +f f cw Y nlle ADF (tom + CEMQU.CEPE) dh 
ctl =f roo f ft 'p(|2 — |) F (2,t, 00 + CP MQ(#), CL P(A) at, 


where k EN,c® = (f',y;), for alli =0,--- ,n and C? = ee Gil; 


an 


3.1.4 Convergence analysis 


In this section, to examine the convergence analysis of our method, we need to recall the 


following useful theorem. 


Theorem 3.2. There exists a constant a > 0 such that, for any function u € H1(I), the 


following estimate holds: 


Pret — ullzaqy S anual acy) 


Proof. For details, see [84] . 


Suppose the following hypotheses: 


4A, A’,B,B’ER,, V¢,7 € [1,1], Vu,u’,v,0' ER, 
|F(¢,7,u,v) -F (¢,7,u',v’)| < Alju-w|+Blv—v', 
JAF (6,7, u,v) — OF (¢,7,u',v’)| < A'Ju—w'| 4+ Blu —v', 
0< V2 (Bll + Blollag) <1. 


Theorem 3.3. Under the assumptions (A), the numerical solution v,(&) converges to the 


exact solution u(€). Furthermore, 


||Un _ Vln) as bn ||o|| (0, 


where 3 is a constant. 
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3.1 Galerkin method 


Proof. Denote Ppv = Un(€) and P,_1v' = v/,(€). Then, from and (3.9), we can write 
vi(€)—v'(6) = 
/ sign( — vie — CF EG, unl) unl) aC + f “al ~ (NEF (E,6, on(Q), v4(6)) dC 
-f sign(e — Qui (KE- DF EG.) ac - f WE = DOF (E,¢, v(6),v'(6)) de 
= J signs — Ong — 6D LF 6, val) e4(6)) a6 ~ FE 6,06), v6) 
+f “le — Cl) [BeF (E,6, on(6), v4 (6)) dC — eF (E,6,0(6), 0G) a6. 


Therefore, under the hypotheses (A), we get 
1 
lon(E) — v'(O) S I IPE — CITA lonlC) — oC) + Blun(¢) — v'(oyIl a6 


1 
a i, IPE — CITA’ un(¢) — v(0)| + B’ lu,.(¢) — v'(Q)I] ac 
< Alp’ llzaq lon — Ull0(1 + Bp’ Ilaaq lle, = U'[ (1 
+ A'llpllz.(y lon = Upc) + Bip llary lon — e'llz~ cy 
Since L® is dense in L?, we achieve 
lee — Olle S V2 (All 'laacy + A'liplliacy) [len — ella 


+ V2 (Bplay + B'lpllxcy) lle — wll aq 


Hence, 


V2 (Alia + A'llallincy 
1= V2 (Blbp'liaw — Biel) 


From the previous Theorem B.J], we conclude that 


lun — U'|ln2q) = llUn — UAlra~M 


ll’m — Ulin < Br“ |olln~, 
where {3 is a constant given by: 


V2 (A lila + A’ lll) 


B=a ! : 
1— V2 (Billa) — Billpllun) 
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3.1 Galerkin method 


n E(n, 2) E(n,3) E(n, 4) E(n, 5) 

3 1.1723e—O1 1.1723e—01 1.1723e—01 1.1719e — 01 
5 2.1376e —03 2.0952e—03 2.0917e—03 2.0915e — 03 
7 5.8607e —04 5.0129e—05 1.6428e—05 1.6271le — 05 
9 5.8485e—04 4.6525e—05 3.2965e—06 8.5842e — 07 


Table 3.1: Numerical results for Example 1. 


3.1.5 Numerical examples 


In order to illustrate the effectiveness of the proposed method, we consider two numerical 
examples for the integro-differential equation (3.1). Let us define the error function as follows: 
E(n, k) = max |v (2) — o» (aI. 


2% 


where u(z) is the exact solution, v,(z) the approximate solution of degree n, z; = —1+ 7s] 
n 


for n = 1---n, and k is the order of the approximation in the algebraic system (3.12). 


Example 1 
Consider the following equation: 
1 
le) = 712) +f /|z —t| cos [z + v(t) +e’ — v'(é)] dé, Vz € [-1,]], 
-1 
v(—1) = -e™, 


with 


f(z) =e — = (V2) + V+ 29) cos(2), 


and the exact solution is v(z) = ze’. 


Example 2 


Consider the following equation: 


v(z) = sin(z) + : (( +2)F+(1- 2)') + I. Jz — t]2./3-+ o(f? + v/(bedt, 


u(—1) =sin(-1), Vz € [-1, 1]. 


such that the exact solution is u(z) = sin(z). 


64 


3.1 Galerkin method 


Exact solution 


* Numerical solution 


25 


a 
Loa) 


Numerical Solution 
a 


o 
a 


-1 -0.5 0 0.5 1 


Figure 3.1: Exact and numerical solutions (Example 1) with n = 2 and k = 2. 


Exact solution 
* Numerical solution 


ois 
on 


Numerical Solution 
—_ 


° 
a 


Figure 3.2: Exact and numerical solutions (Example 1) with n = 3 and k = 2. 
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3.1 Galerkin method 


55'r 


Exact solution 


SP * Numerical solution 


Numerical Solution 


Figure 3.3: The derivative of exact and numerical solutions (Example 1 ) with n = 4 and 


k= 1. 


Exact solution 
* Numerical solution 


Numerical Solution 


-1 a he — —L — 
-1 -0.5 0 0.5 1 
Ss 


Figure 3.4: The derivative of exact and numerical solutions (Example 1) with n = 4 and 


k=3. 
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3.1 Galerkin method 


Exact solution 


08+ * Numerical solution 


Numerical Solution 


Figure 3.5: Exact and numerical solutions (Example 2) with n = 3 and k = 1. 


Numerical Solution 
oO 
T 


-0.2+ 
-0.4/ 
-0.6 + 
-0.8 Exact solution 
* Numerical solution 
-1 
-1 -0.5 0 0.5 1 


Ss 


Figure 3.6: Exact and numerical solutions (Example 2) with n = 4 and k = 1. 
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3.1 Galerkin method 


1A F 


Exact solution 
* Numerical solution 


Numerical Solution 
° ° 
[os] i<e] 
T T 


© 
wy 
T 


0.5 4 4 4 1 4 
0 0.2 0.4 0.6 0.8 1 


Figure 3.7: The derivative of exact and numerical solutions (Example 2) with n = 4 and 


k=l. 


Numerical Solution 


0.55 Exact solution 
* Numerical solution * 
0.5 : ; 
-1 -0.5 0 0.5 1 
Ss 


Figure 3.8: The derivative of exact and numerical solutions (Example 2) with n = 4 and 


k= 2. 
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3.2 Collocation method 


n E(n, 2) E(n,3) E(n, 4) E(n, 5) 

3 3.3368e — 03 3.2190e — 03 3.2189e — 03 3.2188e — 03 
5 5.4483e—04 2.9483e—05 2.6844e—05 2.6746e — 05 
7 5.4679e —04 1.1965e—05 8.6072e—O07 1.3264e — 07 
9 5.4680e—04 1.1966e—05 8.5800e—O07 7.9817e — 08 


Table 3.2: Numerical results for Example 2. 


3.2 Collocation method 


In this section, we present another method to solve the proposed equation, that is the col- 


location method. 


polynomials. 


3.2.1 


Operational integration matrix 


In this method, we seek a numerical solution represented by Laguerre 


The operational matrix of integration for n = 3 is obtained as follows: 


Therefore, 


where 


Similarly, we can find the operational matrix of integration M for an arbitrary value n. 


Yo 
1 
2 


3 


, 
I, 
I, 
[, 


Maxs 


(t)dt 
(t)dt 
(t)dt 
(t)dt 


| 

E 
E 
E 


0 


| S| ER eases 
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C’ P(t)dt = CP MaxsQ(z), 


i 


Beh. 0) 
bes a0 
t,t. 20 
QO “ly 


2-100 0 | Pl), 
-1 100 | P(A), 
0 1 -1 0 | P(z), 


O04 -1| Pe 


3.2 Collocation method 


Let N € N* and consider the following grid points in the interval [—1, 1: 


2 
An = ;=—-l+2 — es 0 he eres A aces 
N {2 +th,h N+1°- 0, ’ ’ \ 


Lemma 3.1. The matrix 


Yo (20) 1 (Zo) Ln (20) 
We = yo(41) 1 (a1) Yn (21) 
Yo (Zn) 1 (Zn) ee) 


is invertible for any n. 
Proof. Let C, be the matrix that contains the coefficients for the Laguerre polynomials of 


degree less than or equals n. We have 


1 
Lis <@: — 
n! 
0 -l ; 
C= : 
_1)r 
ach a) 
n! 
and let V,, be the Vandermonde matrix 
on Bee eee 
by BF eee 
V,= 1 44 1 


It is well known that the matrices C,, and V,, are invertible. We have 


WwW = View 


Hence, the matrix W,, is invertible as well. 


Now, we approximate the function {(z) as 


B(z) = CT P(z), (3:13) 
such that 
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3.2 Collocation method 


BAS Pe), V tS 0pm: 


3.2.2. Method description 


Consider the following integro-differential equation: 


1 
ate) = fle) + J rle—s)¥ (9,81). 8) 4521 yyy 
PHL) = Po: 

By taking the derivatives of (8.14), we get 


: (3.15) 


where 


+1, z>8, 
sign(z—s)=4 0, z=s, 
-l, z<s. 


Let us approximate the unknown function 6’(z) by using (3.13). Then, 


B'(z) = C7 P(z). (3.16) 


By integrating (8.16) with respect to z from 0 to z, we get 


B(z) = Bo + CTMQ(2). (3.17) 
Now, substituting (8.17) and (8.16) into (8.15) yields the following equation: 


1 
C* P(z) = f'(z)+ | sign(z — s)p'(|z — s|)W (z, 8,9 + C’ MP(s),C” P(s)) ds 
ie (3.18) 
: aw 
+ / p(\z — SI) Fy (z, 8, Bo + CT™MP(s), C? P(s)) ds. 
eT 
Collocating the equation (3.18) by the grids points Ay leads to the nonlinear algebraic 


system: 


tL 


3.2 Collocation method 


CTP (z) = f' (%)+ / sign (2; — s) p' (% — 8 |) U (zi, 8, 89 + C” MP(s), C’ P(s)) ds 


1 
: a ; ; (3.19) 
+ f rll) F (5,604 CTMP(s),CTP(S)) ds, 


that can be written as 


OW die ya tO) 2298 (C7), 


where W,, is described in the lemma above and for 7 = 1,...,n, 


di (C") = f' (zi) + / sign (2; — s) p’ (zi — $ |) U (ai, 8, 80 + C’ MP(s), C” P(s)) ds 


1 


L 
+f p(|z — s|) “ (zi, 8, Bo + C7 MP(s),C* P(s)) ds. 
a4 


Therefore, by lemma 3.1] 
=D W,*. 


The nonlinear algebraic system (8.19)is solved by an iterative process. Let us consider 


the following system: 


Cia = [da (C8) ch (CP) ++ sda (CP))] We 
Co = (f' (20) fa) fi Gn), REN. 
By solving the system for a suitable value of k, we can obtain the coefficients C7’. 
Finally, the desired approximation for §’(z) is computed by : 


(3.20) 


3.2.3 Numerical examples 


In order to demonstrate the effectiveness of the method proposed in this section, we provide 
several numerical examples that highlight its accuracy. These examples serve as a way to 
test and validate our method, and compare its results with other existing methods. Through 
these examples, we aim to show that our method can accurately approximate solutions for 


integro-differential equations with a weakly singular kernel. 
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3.2 Collocation method 


E(n, 4) CPUtime 
1.1023e —01 0.045137 
1.553le —02 0.031160 
3.3139e — 03 0.035663 
3.1440e — 04 0.867585 
4.6722e —05 0.045970 


aD oO FF WwW WN] Sd 


Table 3.3: Numerical results (Example 1) 


E(n,3) CPUtime 
4.8956e —03 0.0058 
3.4910e —03 0.1127 
6.7750e —05 ~—-0..1054 
4.8455e —05 ~—(0.1216 
5.4069e — 07 0.1506 


aD oO FF WwW WN] DS 


Table 3.4: Numerical results (Example 2) 


Example 1 


Consider the following equation: 


B(z) = f(z) + I. V/|z — s| cos [z + B(s) +e — B'(s)|ds Vz € [-1,1], 
u(—1) = =e", 


with f(z) = e* — ; (Va —z8+/(1+ 23) cos(z) so that the exact solution is 6(z) = 


ze’. 


Example 2 


Consider the following equation: 


B(z) = sin(z) 4 = ((1 Lz) 4 (1—2)') +f le-sli VS FRCP +R EP as, 
u(—1) =—sin(-1), Vz e€ [-1,]], 


such that the exact solution is G(z) = sin(z). 
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3.2 Collocation method 


Exact solution 


08+ * Numerical solution 


Numerical Solution 


Figure 3.9: Exact and numerical solutions (Example 1), with n = 2 


Numerical Solution 


0.55 Exact solution 
* Numerical solution ‘ 
0.5 
-1 -0.5 0) 0.5 1 


Ss 


Figure 3.10: Exact and numerical solutions (Example 1), with n = 3 
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3.2 Collocation method 


3, 
Exact solution 
* Numerical solution 
25:5 
2 |= 


ai 
ro) 


Numerical Solution 
_ 


0.57 


Figure 3.11: Exact and numerical solutions (Example 2), with n = 2 


Exact solution 
* Numerical solution 


Numerical Solution 


Figure 3.12: Exact and numerical solutions (Example 2), with n = 3 
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3.2 Collocation method 


Discussion 


In this chapter, we have presented two different methods for solving integro-differential equa- 
tions with weakly singular kernels. The efficiency of both methods is demonstrated through 
several examples. Moreover, we can compare advantages and disadvantages of each method. 
The Galerkin method is more accurate but it is quite complex and requires more computa- 
tion. On the other hand, the collocation method is less accurate than the Galerkin method 


but it is simpler and requires less computational complexity. 
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Chapter 4 


Fractional Integro-Differential Equation 


with the Caputo-Fabrizio Sense 


In this chapter, we demonstrate the applicability of our technique in solving fractional integro- 
differential equations with the Caputo-Fabrizio derivative. Specifically, we concentrate on the 


following equation: 


u(z)=9(2) +f K (z,s,u(s), D°u(s)) ds, 
LO) "0; 


(4.1) 


where 0,K,K €C ([0, 1 R’), 0<a < 1,u(z),g9(z) € H'[0, 1], and D® denotes the 


Caputo-Fabrizio derivative of order a. 

We explain the construction of the fractional operational matrix of integration utilizing 
Hermite wavelets. We then apply the collocation technique followed by the iterative method 
to approximate the numerical solution. To verify the accuracy of our process, we present 


several computational examples. 
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4.1 Hermite wavelets 


Figure 4.1: Hermite wavelets for n = 0,1, 2,3. 
4.1 Hermite wavelets 


The Hermite wavelets are given as follows : 


k+1 ‘ : 
22 1~—1 2 


9k-1 Zz 9k-1? 


AAs 91), 
i= A ee ) 


a 


otherwise , 


where i = 1,2,...,2*-1.k > 0 is an integer number, 7 = 0,1,2,...,n—1, H, is the Hermite 


polynomial of degree 7. 


Any function u(z) in L2,(IR) can be written in the following form: 


u(z) = Ss" > cig Vigl2), 


i=1 j=0 


where c;,; = (u, v;,;) with (.,.) being the scalar product in the Hilbert space L7( 
we get the approximate function for u(z) by truncating the series (4.2) as follows: 


Qk-1 n-1 


ile) = S- », Ciittig(z) = C7? P(z), 


i=1 j=0 


where P(z) and C” are 2°~'n x 1 matrices: 
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(4.2) 


R). Then, 


(4.3) 


4.2 Fractional operational integration matrix 


T 
P(z) = [w10, W115 tee Wie ais 20, W2,1, sae i a see »Wok-1.0, Wor-14 siece , Yok-1,n—1| ) 


and 


T 
C = [e10. C115 are »C1n—-1) C105 C215 sang »C2n—1) cana »Cok-1 9, arta , Cok-1 1 | a 


4.2 Fractional operational integration matrix 


If k = 1, then both P(z) and C” would be: 


Ge = [a0, Olfisncces Ona! 5 
PE) = [o(z), ea(z), ae Waal?) : 


W,, denotes a matrix comprising the coefficients associated with the Hermite wavelets: 


2-4 4 = --- Hy 1(-1) 
0 8 32 : 
1 
W, nl oie : 5 
Th O90 332 

0 0 O- =: See 

and 
Fila) eee ae” YN, ek Pale al) eae) 

We have: 


Pal Z) = Znl2)W x 


Consider an integral matrix denoted by N within the canonical polynomial basis: 


Oc 0p i. ee < 30 
OG Dae GO! see - 210 
L= | 0°. 0. Te se <0 
Ort JO 0 he fe 


Subsequently, the operational integration matrix in the Hermite wavelets basis is given as 


follows: 


| CTP, (y)dy = CT MPraa(2), 
0 
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4.3 Method description 


where 


M =WLWast. 


Furthermore, employing the previously defined notation in (1.1), we can express the 


fractional operational integration matrix in the following manner: 


L* (CP, (z)) = C*[aM + (1- a) F)Pasi(z) = C7Qn(z), (4.4) 
10-5 0 
0 


where Q),(z) = [(l-a)F +aM|P,4i1(z) and F= 


4.3. Method description 


Consider the following Fredholm integro-differential equation: 


ulz) = fl2) +f K (z,8,u(s), D°u(s)) dy, u(0) = 0. (4.5) 


By taking the derivative of both sides of equation (4.5) and using the Caputo-Fabrizio 


derivative of order a, we obtain 
1 
Dus) =D" (2) +f DEK (z, 8, u(s), D°u(s)) ds. (4.6) 
0 
We approach the unknown function D°u(z) by using the formula (4.3) 


D°ulz) = CP, (2). (4.7) 


To approximate the unknown function u(z), we perform an integration of (4.7) with the 


help of the operational matrix of fractional integration mentioned earlier (4.4). This process 


yields: 
Un(Z) = I* (CTP, (z 
(z) ( (z)) (4.8) 
=O Oz): 
Now, substitute and into to obtain: 
1 
GTP, (2) =D" q{2) +f Dek (24,0 O,6).C" Ps) ) de (4.9) 
0 
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4.4 Numerical experiments 


2a 1 
2(n + 1) 
equation (4.9). This results in the formation of the subsequent nonlinear algebraic system: 


Using the grid points z; = , where i € N andi < n—1, we apply collocation to 


CTA=D", (4.10) 


such that 
D® = [do, di, +++ ,dn—1] 


where 
1 
d; = D%g (z%) + | Dok (2:,8,C7Qn(s),C* Pa(s)) ds, fori =0,--- ,n—1, 
0 


and 


Wo (20) o(z1) +++ Po (Zn-1) 
Wi(zo) = Wi (41) +++ U1 (n-1) 


Yn—1 (Zo) Pn-1(%) +++ Yn-1 (%n-1) 
We apply the iterative method to solve the system (4.10). To do this, we introduce the 
following system: 
Cat =D. 
where 
Dy = dh ar ade | 5 


and 


1 
dk =D%q (x) + i DEK (zi, 8, Ck Qn(s), Ck Pa(s)) ds. 
0 


For a suitable value of k, we find the vector C.. then substitute the coefficients of Gr 
into (4.8) to compute the approximate solution to equation (4.5). 
4.4 Numerical experiments 


Here, some illustrative experiments are included to demonstrate the efficiency of our method. 


We introduce the error as follows: 


En, = max |un(z) — u(z%)], 
i=0,n—1 


where u(z) represents the true solution, u,,(z) stands for the approximate solution, n signifies 


the degree of Hermite wavelets. 
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4.4 Numerical experiments 


First experiment 


Let be given the fractional Fredholm integro-differential equation: 


Vz € [0,1], u(z) = g(z) +f In E (cos(s) — ze~**) + =u(s) — D°*u(s)} ds, 


where 


ieee Sa r 2] 7 


3 
5" 
The exact solution to this equation takes the form of u(z) = sin(z), when the fractional 


order of differentiation is a = 0.5. 


Second experiment 


Consider the following equation: 


7 ; sin(z + s) 
te Se) i; 1+ D°u(s) + ey mee Al 


In this context, we have g(z) = e * — 1+ cos(1 + z) — cos(z), and the exact solution is 


expressed as u(z) =e * — 1 if the fractional order of differentiation is a = 0.75. 


Third experiment 


We have the following equation: 


5 1 
Vz € [0,1], u(z) = =5f + | 27/12 — 3e-8 + 2D2u(s) + u(s)ds, 
0 


2 
such that the order of derivation is a = 3? and the exact solution of the equation is u(z) = 27. 


Fourth experiment 


Let the following equation: 


z t 2 A s s 
Vz € [0,1], u(z) = ze 3 — cos(z) + | sin (2 —s+ 5° — 3% + e3u(s) + ciDu(s)) ds, 
0 
(0) =O: 


1 2 
The order for this example is a = mt and the exact solution is u(z) = ze~ 3 — cos(z). 
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4.4 Numerical experiments 


n 3 4 5 6 7 
n | 6.988E — 05 | 1.433E — 06 | 4.934E — 08 | 6.723e — 10 | 5.444E — 10 
Table 4.1: Numerical results (First experiment). 
n 3 zi 5 6 vi 
Ey, | 7.949 — 04 | 2.759E — 05 | 7.617E — 07 | 1.741E — 08 | 8.503E — 09 
Table 4.2: Numerical results (Second experiment) 
n 3 4 5 6 7 
EL, | 4.402E — 04 | 5.399E — 05 | 5.730E — 06 | 4.661E — 07 | 4.355E — 08 
Table 4.3: Numerical results (Third experiment) . 
n 3 4 5 6 c 
Ey, | 1.301E — 05 | 1.762E — 07 | 1.065E — 08 | 7.066E — 11 | 4.737E — 11 


Figure 4.2: Exact and approximate solutions (First experiment), n = 7. 


Table 4.4: Numerical results (Fourth experiment). 


Numerical Solution 


-0.7 


Exact solution 
* Numerical solution 


0.2 0.4 
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0.6 0.8 


4.4 Numerical experiments 


Exact solution 
0.8; * — Numerical solution 


Numerical Solution 


Figure 4.3: Exact and approximate solutions (Second experiment), n = 7. 


Exact solution 
* Numerical solution 


Numerical Solution 


Figure 4.4: Exact and approximate solutions (Third experiment) , n = 7. 
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4.4 Numerical experiments 


0.85 


Exact solution 


* Numerical solution 
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T 


Numerical Solution 
Oo Oo 
w bh 


© 
i) 
1 


0.1 F 


Figure 4.5: Exact and approximate solutions (Fourth experiment) n = 7. 


Results interpretation 


Tables and [£4 display the error E,, for various degrees n, illustrating that 
the method’s performance improves as n increases. Additionally, Figures 4 and 
present graphs depicting both the exact and approximate solutions, showcasing their 
remarkable similarity. Thus, the aforementioned examples serve as compelling evidence of 


the efficiency and validity of our numerical approach. 
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Conclusion 


This thesis has been focused on an in-depth investigation of various types of integrodifferen- 
tial equations. Specifically, we have studied equations with regular kernels, weakly singular 
kernels, and those in the fractional case. To solve these equations, we have employed the 
popular projection method along with classical orthogonal polynomials. By applying this ap- 
proach, we are able to transform in each case the main equations into a nonlinear algebraic 
system, that can then be solved using the Picard successive approximations. To validate 
the effectiveness and accuracy of our proposed methods, we have also developed algorithms 
using the Matlab platform. This allowed us to implement our methods and present numerical 
examples that demonstrate their applicability and performance. Through our comprehensive 
analysis and numerical examples, we have shown that our proposed approach can accurately 


and efficiently solve a wide range of integro-differential equations. 
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